Conceptual analysis of the combined effects of vaccination, therapeutic actions, and human subjection to physical constraint in reducing the prevalence of COVID-19 using the homotopy perturbation method

Background The COVID-19 pandemic has put the world's survival in jeopardy. Although the virus has been contained in certain parts of the world after causing so much grief, the risk of it emerging in the future should not be overlooked because its existence cannot be shown to be completely eradicated. Results This study investigates the impact of vaccination, therapeutic actions, and compliance rate of individuals to physical limitations in a newly developed SEIQR mathematical model of COVID-19. A qualitative investigation was conducted on the mathematical model, which included validating its positivity, existence, uniqueness, and boundedness. The disease-free and endemic equilibria were found, and the basic reproduction number was derived and utilized to examine the mathematical model's local and global stability. The mathematical model's sensitivity index was calculated equally, and the homotopy perturbation method was utilized to derive the estimated result of each compartment of the model. Numerical simulation carried out using Maple 18 software reveals that the COVID-19 virus's prevalence might be lowered if the actions proposed in this study are applied. Conclusion It is the collective responsibility of all individuals to fight for the survival of the human race against COVID-19. We urged that all persons, including the government, researchers, and health-care personnel, use the findings of this research to remove the presence of the dangerous COVID-19 virus.


Background
COVID-19 infection has spread to every continent. New coronavirus hotspots were detected in Wuhan, China, in December 2019. When numerous people were admitted to the hospital in late December 2019, it was clear that the pandemic had begun. By mid-July 2020, the virus had infected over 213 countries, causing 15,969,465 infections and 643,390 fatalities [1]. The WHO discovered that the virus may be breathed in through normal breathing, resulting in new infections. COVID-19 has a 2 to 14 days of incubation period, with about 97.5 percent of infected individuals presenting symptoms 11 days after infection [2][3][4].
Several researchers and scientists have published significant and successful research works to aid in the eradication and control of COVID-19. Among these is the study reported in [5], in which the fear impact of media is divided into two categories, namely the fear of infection toll and the fear of death toll, and investigated using a proposed mathematical model. Additionally, the influence of the convex incidence rate on multiple COVID-19 transmissions was explored in [6]. According to their findings, the double exposure of vulnerable individuals can result in projected occurrences of COVID-19 infection in the general population. A novel assessment of the COVID-19 transmission dynamics was reported in [7]. The model looked at the effects of the convex incidence rate. To stabilize the pace of disease reduction in their system, a variety of strategies were utilized. Eventually, the dynamics of the disease are greatly impacted as huge drops in infection rates are observed when the rate of immigration and population mixing in the afflicted area are taken into account.
Vaccination, therapy, and human compliance to restriction of physical interactions are crucial in avoiding illness development associated with COVID-19. Several researchers have developed mathematical models to solely evaluate the impact of each of these factors in eradicating the spread of COVID-19. The results of a twostage COVID-19 vaccination research were published in [8]. According to their findings, well-implemented immunization programs can prevent the spread of COVID-19. In [9], research was done on the SEIRS model for COVID-19 that captured saturation incidence and treatment response for the global study. In their investigation, they came to the conclusion that, even if treatment can reduce the spread of COVID-19, lowering the effective contact rate is the best way for people and the government to do so. Additionally, a mathematical model that can evaluate the effect of mask use on the general populace was proposed in [10], and a SIR COVID-19 pandemic model was investigated using statistical methods in [11].
The course of diseases may usually be accurately and realistically predicted by performing a numerical simulation of epidemic models. To forecast the COVID-19 virus's future existence, a mathematical model of the virus was built, investigated, and simulated in [12] using actual data from Pakistan. Researchers in [13] looked at the logistic growth model that might be used to estimate the scope of the COVID-19 epidemic. In [14], a mathematical analysis of a stochastic model for coronavirus propagation was reported. A collocation approach based on Legendre polynomials was used to achieve the numerical solution of this system, and simulations were used to review the results and pandemic model findings, and a recent research presented in [15] features the simulation of fractional-order Caputo's derivative on a coronavirus disease model using the Laplace-Adomian decomposition method.
Since these mathematical models are typically nonlinear, strong numerical techniques, such as the homotopy perturbation approach offered by [16], are often required to compute their approximate solution. This homotopy perturbation approach was used by researchers in [17] to solve an epidemic model of EIAV infection. Their findings demonstrate how effective the method is at resolving coupled nonlinear differential equations. This same method has also been applied by [18], where the effect of the disease transmission coefficient on a disease-induced death seizure epidemic model was simulated.
As of date, no research has yet been published that addresses how to restrict the spread of COVID-19 while simultaneously taking therapeutic interventions, immunization, and human cooperation to limit physical interaction into account. The reasons for including these factors in our research are discussed in [8,19,20]. Using extracted COVID-19 data in the literature and real-life data obtained from the Nigeria Centre for Disease Control (NCDC) [21], we shall conduct theoretical and numerical analysis on a newly proposed mathematical model of COVID-19 transmission dynamics, which is a modification of the mathematical model in [22], where the dynamics of COVID-19 in relation to isolation class was studied.

Model design
The predicted dynamics of COVID-19 transmission in the population are shown in the schematic diagram below.
This diagram might be read and expressed mathematically such that

Descriptions
Equation (1) represents a deterministic mathematical model of COVID-19; the population of the model is divided into five classes, namely the susceptible class S(t) , exposed class E(t) , infected class I(t) , quarantine class Q(t) and recovered class R(t) . The therapeutic response action T for the sick person I is a piecewise linear function defined in [19]. It varies on the interval 0 ≤ r < 1 and is related to the prevalence of infected individuals. Another control parameter under study is "c. " It combines curfews, social distancing, and other physical contact restriction measures which are examined on the model to possibly limit the rise of the effective contact rate between the susceptible and infected populations. The immunization of vulnerable persons is represented by parameter ρ . We theoretically assume that all susceptible people, including the current and newly recruited members of the class, are administered the vaccine in order to assess its influence on the prevalence of the disease.
For biological reasons, physical contact between the classes of vulnerable and sick individuals should be limited to reduce the disease transmission coefficient. Therefore, it makes sense to create a function f (β, c) in such a way that:

As such
(i) f (β, 0) = β , i.e., without physical restraint, interaction between the classes will continue and the viral coefficient will spread progressively and infection will persist in the system. (ii) f (β, 1) = 0 , i.e., the model's classes won't interact and disease transmission in the system will stop.
Thus, if there is no infected individual in the system, there is no requirement to prohibit physical contact among the populace. (1) Since the therapeutic action T on infected individual I(t) I is rI(t) , the system of Eq. (1) becomes: Other factors involved in the dynamics of the disease transmission are described in Tables 1 and 2.

Qualitative analysis 2.3.1 Existence and uniqueness
A Lipchitz criterion will be employed to ensure that the solution exists and is unique. Thus from Eq. (2), let: Following the criterion, we obtain the system's partial derivatives.

Nonnegativity of invariant region
The invariant area establishes the viable region for the solution of an epidemiological model, often known as the limit or range for the solution of the model. Any solution outside of this range is neither Since the human population varies throughout time, which yields

Nonnegativity of Solution
Separating the variables, Integrating both sides and applying the initial conditions, Following the same procedure, we demonstrate the positivity of the other classes. Thus, The solutions are positive, and this completes the proof.
After satisfying all of the fundamental requirements for an epidemiology model, we conclude that the suggested model is appropriate for studying the dynamics of COVID-19 in the general population.

Analysis of equilibrium states 2.4.1 Disease-free equilibrium (DFE)
By setting the right side of (2) equal to zero, we may determine the equilibrium point devoid of sickness, i.e., We get the following results by solving the equations that emerge from substituting the disease-related classes E, I, and Q to zero: Equation (9) represents the disease-free states of the proposed model.

Endemic equilibrium state
Endemic equilibrium, also known as a nonzero equilibrium condition, occurs when a disease persists in a population. Contrary to the disease-free state, For brevity of terms, we let: Thus, the endemic equilibrium state yields

Basic reproduction number
As defined in [23], the estimated number of secondary cases that can be directly caused by a single case in a population that has reached the peak of an illness is known as the basic reproduction number R 0 .
To compute the basic reproduction number ( R 0 ), we take into account the following two disease-manifested compartments.

Local stability analysis of disease-free equilibrium
While investigating the stability of the model's equilibria, the following results are demonstrated: Proof Consider the Jacobian of Eq. (2), which is given by: The eigenvalues are then calculated, and the results are as follows: We infer that the disease-free equilibrium exists since the eigenvalues are all negative.

Global stability at disease-free equilibrium
Lemma 2 Applying the Lyapunov function approach to the global stability of the disease-free equilibrium state of the model, . According to LaSalle's invariance principle, the disease-free equilibrium E 0 is globally asymptotically stable if R 0 < 1.

Sensitivity Analysis
We used the sensitivity index analysis to assess the resilience of the model parameter based on the fundamental reproduction number. The sensitivity indices of each parameter can be computed with the following formula: According to the findings, the three control variables yield low sensitivity indices. They are negative, but given how close to unity they are, it is hard to underestimate their effectiveness in lowering the basic reproduction number (Table 3).

Homotopy perturbation method
One of the aims of this research is to do numerical simulations on the mathematical model. To do so, we must first establish an approximation solution to the mathematical model. In pursuit of this, we explore the methods of applying the homotopy perturbation. The method's analysis is explained more below. Hence, consider the following differential equation: Subject to the boundary condition Operator denotes the differential operator, the boundary operator is , k(r) is an analytic function, the boundary of the domain is denoted by , and ω n is the normal vector derivative drawn externally from . We can split the operator �(ω) into two parts such that    We can naturally assume the solution (28) as a power series such that Such that evaluating (30) using (33), and comparing coefficients of equal powers of p.
The values of f 0 (t), f 1 (t), f 2 (t) are obtained by solving the resulting ordinary differential equations. Thus, the approximate solution of (25) is:

Numerical simulation
In this part, we use the homotopy perturbation approach to conduct the numerical simulation that produces the SEIQR epidemic model's approximate solution. Constructing a homotopy for (1), The approximate solution of (1) can be assumed as: Substituting (36)   S(t) = s 0 (t) + ps 1 (t) + p 2 s 2 (t) + . . . p n s n (t) E(t) = e 0 (t) + pe 1 (t) + p 2 e 2 (t) + . . . p n s n (t) The coefficients of p 2 equally yield: The second approximations are obtained by solving these equations.
This continues till the desired number of iteration is computed. For the purpose of this research, we performed three iterations. Due to cumbersome number of terms, we present the iteration code instead of the computed results.

Results
The solution for each class is obtained by taking the sum of its obtained approximations such that:

Convergence of HPM
The convergence of the homotopy perturbation approach strictly depends on the contraction of the approximation solution to the exact solution [25].

Theorem 3
Let there exist a mapping k : m → n defined on two Banach spaces m , n for all s, t ∈ m then k(s) − k(t) n ≤ δ�s − t� m , 0 < δ < 1 such that the +· · · sequence s r+1 = k n (s 0 ) = τ (s 0 ) for some s 0 ∈ m which converges to a unique fixed point k [25].
Proof We consider a Picard sequence s r+1 = k(s r ) ⊆ n to prove the theorem. It is required to show that s r is convergent in n for all r ≥ ν such that �s r − s v � ≤ �s r − s r+1 � +�s r+1 − s r+2 � + +�s r+2 − s r+3 � + · · · +�s r−1 − s v � . The proof is defined by applying mathematical induction on the contractive property of contraction c such that This proves that (s r ) is convergent in n and through completeness of n , we can find ω ∈ n : lim r→∞ (s r ) = ω ∈ n . Clearly, contraction C ensures the continuity of k . Thus, ω = lim r→∞ s r+1 = s v .

Lemma 3
The proposed mathematical model's convergence of s r to s v cannot be asserted since no exact solution to the model exists. Table 4 Third-order maximum truncation error for t ∈ [0, 1]  Table 5 Error range of each class at t ∈ [0, 1] Proof The rate of convergence of s r depends on the following conditions as suggested by He [16,26]: (i) The nth-order derivative of f (t) with respect to t in (33) must be small as parameter p may be relatively large as p → 1 (ii) ii. ℓ −1 ∂f (t) ∂v < 1 so that the series converge.
This rate is examined in the following section.

Error analysis
In numerical computations of series solutions, error analysis is crucial. In this section, we utilize the convergence requirements from Lemma 1 to perform an error analysis on the approximate results of each compartment of the model using Maclaurin's estimation of maximum truncation error at predefined baseline values of the parameters given in Table 2.
Definition The truncation error of an nth-order Maclaurin's polynomial of a function f (t) with n + 1 derivatives on the interval |t| ≤ r and in the interval f n+1 (t) ≤ M has the following bound as maximum error E n (t) = max To compute the third-order maximum truncation error of the model's class series results on the interval [0, 1] , we construct the following functions: The error analysis of the approximate series results demonstrates that the model is validated for actual data simulation since the error generated by the third-order approximation is negligible. We observed that the maximum error is recorded in the susceptible and exposed compartments. This happens because the two classes contain nonlinear variables (Tables 4 and  5; Fig. 1).

Numerical simulation
In this section, we execute numerical simulations of the mathematical model to rigorously investigate the influence of the incorporated parameters ρ, c &r on the interval [0, 1) to ensure the effectiveness of the regulatory factors integrated in the research. Figures 2,3,4,5,6,7,8,9, and 10 depict the simulation process's outcomes graphically. Figure 12 reveals the various responses of the exposed and infected human population to combine the influence of the incorporated control parameters.

Validity of the model
In this section, we assess the model's validity by numerically projecting active infected cases for 30 days using the baseline parameters given in Table 2 and comparing the outcome with a time series plot of daily data of COVID-19 reported by the Nigeria Centre for Disease Control between December 1 and 31, 2020 [21]. The model's validity is evidenced by the fact that it achieves its intended objective. The validation result of the proposed model in Fig. 11 demonstrates that the proliferation in active infected cases when the study's incorporated measures are not used correlates with the time series plot of real COVID-19 data.

Discussion
The results of the numerical simulations are discussed in this section. The restriction on population-wide physical contact, extensive imposition of immunization, and the intense care given to the sick patients were all clearly examined and recorded to have increased the number of susceptible and recovered individuals in a COVID-19-enraged population. Remarkably, we observed that the findings of this analysis are consistent with those of earlier research [9], which focused on the impact of piecewise linear therapy on a novel SEIRS model of coronavirus (COVID- 19) illness. The maximal imposition of therapy is also included in our study, and it was demonstrated that the count of recovered individuals peaked when therapy is imposed at a maximum level. The simulation results shown in Figs. 5 and 6 revealed that the therapy parameter incorporated in our model exhibits similar responses in recovered and infected patients as obtained earlier in [9]. Hence, therapy is crucial in limiting the transmission of the COVID-19, and we propose that therapy-associated factors such as treatment and medications be effectively levied on infected people in order to increase the number of recovered COVID-19 patients. Figures 2, 3, and 4 examine the impact of physical limitation compliance "c" on the SEIR class. In a similar but distinct study reported in [20], raising the probability of an individual's positive reaction to the curfew and social distancing lowers the basic reproduction number below one, and this will slow the progression of the disease. there is no infection transmission between the two groups and no progression from susceptible to exposed and infected, indicating that system (1) will eventually approach asymptotically stable disease-free point. Figures 3 and 4 assert that restricting physical contact would substantially lower the number of exposed and sick persons; hence, if the masses failed to comply with the rule of curfew, social distancing, etc., system (1)'s trajectory may approach the endemic point, and all individuals in the recovered class will be wiped out. Although increasing the response rate of an individual's body contact in a disease-manifested system cannot totally prevent the disease's spread, it can reduce it to a bare minimum; hence, we explore the effect of vaccination in reducing infectious disease transmission. Figures 7, 8, 9, and 10 depict the response rates of susceptible, exposed, infected, and recovered population systems to immunization.
As revealed in an earlier study [8], immunizing vulnerable individuals minimizes the number of exposed and sick people. Vaccination also showed significant effects in limiting infection propagation in the current study, as recovered individuals were observed to peak when vaccination was fully provided. As a consequence, the disease will be eradicated by the proposed system. The collective impact of the three control mechanisms on the mathematical model is depicted in Fig. 12. In the absence of the three control factors in Fig. 12, it was observed that active cases of exposure and infections of COVID-19 keep projecting. Contrarily, this increasing population of the two classes was observed to be drastically lowered to minimum when the control mechanisms are utilized. Also, the homotopy perturbation approach used in this work has been used to analytically solve and mathematically explore a suggested smoking model utilizing the Caputo-Fabrizio operator in [27]. The approach also effectively generated the solution of our proposed mathematical model since it yielded a stable result which unconditionally describes the dynamics of coronavirus disease, as seen in our simulation results. If time is increased, as seen in several of the graphs, negative population figures will be recorded. An identical dynamic was identified in [28], when the Laplace-Adomian decomposition approach was used to analyze the effect of fractional-order on a smoking model. Because of the small data sets employed, a short time interval was used to prevent a negative population of classes, which would be unrealistic; thus, we propose that for a longer time interval, big values of the initial data should be used so that the concerned population remains positive. Furthermore, since there are several factors associated with the dynamics of disease transmission during pandemic, short time intervals should be considered for authentic results due to different sensitivities of the model parameters to time.

Conclusions
Epidemic infections, such as COVID-19, cannot be totally eradicated by a single factor. The objectives of this research were achieved as the numerical simulation analysis unconditionally showed the dynamics of COVID-19 spread under different control situations. We also showed that the combined application of vaccine, therapy, and human compliance factors to COVID-19 curfew and social distancing rules has an appreciable effect in stopping the disease's global existence. The mathematical model's validity and error analysis showed that the model could be applied for further studies on coronavirus disease as it depicts the behavior of real-life COVID-19 situations.
To eradicate the COVID-19 virus, we recommend that governments around the world give immunizations to their vulnerable populations on high range. Furthermore, hospitals should provide superior health-care facilities as well as advanced life-saving technology capable of treating sick people flawlessly. Furthermore, it is critical to underline that governments' efforts to restrict the transmission of the virus can only yield beneficial results if the populace has a receptive attitude toward avoiding disease contact. To constantly eliminate the disease's existence, it is thus recommended that everyone obey curfew, isolation, and the use of nose masks, among other things.

Recommendations
Although this study demonstrated the diverse reactions of an endangered population to therapy, physical contact limitation, and immunization in a COVID-19-exhibited environment, future research could incorporate introduction of second-stage immunization strategy as well as incidence rate capturing fear of disease transmission. Also, fractional-order analysis of the current mathematical model employing fractional operators like Caputo [29], Caputo-Fabrizio [30], or Atangana-Baleanu [31] might also be useful in examining the disease's dynamical trend on real data.